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I. 


Introduction 


To  date,  the  energetic  plasmas  responsible  for  che  terrestrial  aurora 
have  been  probed  directly  with  in  situ  spectral  measurements  of  the  precipitat¬ 
ing  electron  flux,  while  the  morphology  and  temporal  behavior  of  the  aurora 
have  been  monitored  with  ground-based  and  space-based  imagery.  Space-based  X- 
ray  imaging  offers  a  means  for  combining  spectral  and  spatial  measurements  of 
the  aurora  by  detecting  X-ray  fluxes  which  arise  from  bremsstrahlung  emission 
by  the  precipitating  auroral  electrons.  Of  course,  the  bremsstrahlung  X-ray 
spectrum  must  be  interpreted  properly  in  order  to  infer  the  spectral  parame¬ 
ters  of  the  precipitating  electrons.  Quantitative  interpretation  of  brems¬ 
strahlung  X-ray  spectra  is  possible  since  the  production  and  transport  of  >  1 
keV  bremsstrahlung  X-rays  in  the  atmosphere  are  reasonably  well  understood. 
Also,  observations  of  bremsstrahlung  X-rays  have  been  accomplished  (e.g., 
Mizera  et  al. ,  1978;  Imhof  et  al. ,  1974,  1975a,  b),  and  the  most  recent  obser¬ 
vations  have  provided  two-dimensional  images  of  the  X-ray  flux  along  with 
spectral  measurements  (e.g.,  Mizera  et  al. ,  1984;  Imhof  et  al.  ,  1982). 

In  practice,  the  use  of  bremsstrahlung  spectral  measurements  to  infer 
incident  electron  spectral  parameters  is  hampered  by  instrumental  limitations 
and  low  signal-to-noise  ratio.  Thus,  a  spectral  deconvolution  scheme  which 
makes  the  most  effective  use  of  available  data  is  required.  The  maximum 
entropy  formalism  (Jaynes,  1957)  is  a  statistical  technique  which  provides  a 
statistically  optimal  estimator  of  a  continuous  quantity  (in  this  case,  the 
incident  electron  spectrum)  when  only  a  discrete  set  of  integral  functions  of 
that  quantity  are  known  (namely,  the  observed  X-ray  spectrum).  Maximum  entro¬ 
py  analysis  has  been  applied  previously  to  image  restoration  problems  (e.g., 
Andrews  and  Hunt,  1977;  Gull  and  Daniell,  1978,  Wilczek  and  Drapatz,  1983) 
where,  for  example,  brightness  Images  of  astronomical  objects  appear  noisy  and 


blurred  due  to  inadequate  detection  equipment.  The  bremsstrahlung  spectral 
deconvolution  problem  can  be  compared  to  noisy  and  blurred  image  restoration, 
since  the  electron-atmospheric  interaction  effectively  causes  "blurring"  of 
the  original  electron  spectral  signature. 

A  numerical  maximum-entropy  deconvolution  scheme  has  been  developed  which 
is  specific  to  the  bremsstrahlung  problem.  The  technique  includes  a  Chi- 
squared  limitation  on  errors  arising  from  statistical  fluctuations  in  the  X- 
ray  flux  measurements.  The  numerical  technique  has  been  applied  to  data 
acquired  by  the  Aerospace  X-ray  spectrometer  on  the  polar-orbiting  DMSP-F6 
spacecraft. 

Section  II  of  this  report  describes  the  fornxilation  of  the  appropriate 
bremsstrahlung  production  function  (i.e.,  the  "blurring"  or  "response”  func¬ 
tion).  Section  III  describes  the  maximum  entropy  constraints  and  derives  a 
linear  system  of  equations  which  comprise  an  interactive  solution  to  the 
problem.  Section  IV  applies  the  technique  to  some  ideal  and  actual  examples, 
and  Section  V  summarizes  the  results. 


II.  Bremsstrahlung  Production 


It  is  necessary  to  determine  an  appropriate  bremsstrahlung  production 
function  $(X,  T,  n)  which  approximately  satisfies  the  following  expression  for 
the  observed  X-ray  spectrum: 

T 

fx(X,  n)  »  /  MAX  dT  fe(T)  •  <KX,  T,  n)  +  e(X)  (1) 

x 

where  the  integral  for  the  observed  flux  fx  as  a  function  of  X-ray  energy  X 
and  aspect  angle  n  (see  Figure  1)  must  be  carried  to  an  incident  electron 
energy  TMAX  such  that  the  contribution  to  observed  X-ray  flux  from  incident 
electron  flux  feo(TMAX^  is  neg11-?1-1316*  The  factor  e(X)  represents  errors  or 
variations  which  arise  from  processes  not  described  by  the  production  function 
Later  we  assume  that  the  errors  e(X)  are  not  biased  and  represent,  for 
example,  fluctuations  due  tu  inadequate  counting  statistics.  In  practice  the 
production  integral  (1)  is  replaced  with  a  summation  over  energy  levels 
which  are  selected  for  evaluation.  That  Is 

n 

f  ,  ( n)  *  E  f  .•$..(  n)  ♦  e  (2) 

xk  .  ,  e  i  i , k  k 

l*  1  * 

The  bremsstrahlung  production  kernel  $  can  be  estimated  based  on 

I ,  k. 

analytical  treatment  of  bremsstrahlung  production  and  transport  in  the  atmo¬ 
sphere  (e.g.,  Luhmann,  1976;  1977)  or  from  the  results  of  Monte  Carlo  calcula¬ 
tions  (e.g.,  Seltzer  and  Berger,  1974).  We  review  here  an  analytical  tech¬ 
nique  for  estimating  observable  X-ray  flux  at  satellite  altitudes  in  order  to 
demonstrate  some  important  physical  considerations  (limb  brightening,  for 
example).  In  any  case,  any  appropriately  derived  bremsstrahlung  production 

function  $  can  be  substituted  into  the  maximum  entropy  formulation  de- 

I ,  k 

scribed  in  Section  III. 
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The  bremsstrahlung  production  function  <t»  represents  the  integrated  X- 

ray  productivity  of  an  incident  electron  beam  at  energy  over  an  emission 

range  0  <  z  <  z„  where  2  is  the  atmospheric  depth  (gra/cm^).  The  form  of  & 

c  i ,  k 

must  incorporate  both  electron  energy  degradation  over  the  entire  emission 
range  and  X-ray  absorption  over  the  escape  path  length  z  sec  n,  where  the 
angle  n  is  the  emission  aspect  angle  with  respect  to  the  zenith.  A  plane 
parallel  geometry  can  be  assumed,  as  shown  in  Figure  1.  Noce  that  f  (T)  must 
be  assumed  uniform  over  the  "sensed"  horizontal  area. 


Luhmann  (1976)  treated  the  spatial  evolution  of  the  energy  spectrum  of 
incident  electrons  on  the  atmosphere  by  showing  that  a  solution  to  the  Fokker- 
Planck  equation  for  the  electron  distribution  functions  has  the  form 


fg(z ,E) 


2*  dT/dS(T) 
v  dT/dS(E)  eou; 


(3) 


where 


Y-  3/4 

dl7dS(Y)  *  7?TT 


and 


T  - 


r! 
1 B 


(z+BE 


7/4 


(4) 


(5) 


for  a  beam.  Here  Y  is  a  dummy  variable,  E  is  the  energy  of  the  degraded  beam 
and  B  *  4.57  x  The  range  of  the  beam  is  thus  zQ  *  BT'14,  however  it  is 

more  important  to  consider  the  emission  range  of  the  beam,  beyond  which  an 
electron  beam  with  InitiaL  energy  T  does  not  contribute  to  brerass t rah  lung 
production  at  photon  energy  X.  This  range  is  such  that  E  *  X  in  Eq .  (5),  or 


1 1 


B  (V  /4  -  X7/4  ) 


(6) 


Similarly,  one  can  determine  the  minimum  incident  electron  energy  which  can 
contribute  to  bremsstrahlung  emission  at  energy  X  for  atmospheric  depth  z, 
which  is 


I  -  [i  (.  ♦  bx7/4)] 


The  flux  spectrum  of  electrons  at  depth  0  <  z  <  z „  is  then  just 


fe  (z,E)  -  (i)  feo(T)du 


where  Au  is  the  cosine  of  the  beam  width  in  pitch  angle. 

The  bremsstrahlung  production  rate  due  to  an  electron  of  energy  E  is 
determined  by  the  Bethe-Heitler  differential  cross  section  (Heitler,  1954), 
which  for  an  air  target  of  mean  atomic  weight  M  takes  the  form 


1  d  °B 
M  dX 


fv  .x  1.75  ,  r(  «He  +  /~E=X)Z  1 
(X,E)  --ysr-  in  I - - -  j 


with  units  cm^/g  keV.  This  cross  section  is  plotted  in  Figure  2  for  selected 
electron  energies.  Note  that  high  energy  electrons  contribute  significantly 
to  X-ray  production  at  much  lower  energies. 

The  bremsstrahlung  production  rate  at  atmospheric  depth  z  is  thus 


n  ,  «  d  a„  c  3/4 

V2’"  "  V  I  dT—  lx'E)  f.o<T) 

c 


where  the  factor  of  Q  arises  from  an  integration  over  the  angular  distribution 

of  bremsstrahlung  emission  with  respect  to  the  direction  of  the  electron 

beam.  Since  we  are  concerned  with  auroral  energies  only,  the  slight  energy 

dependence  in  the  emission  angular  distribution  (Rees,  1964)  has  been  Ignored 

2 

here,  and  for  an  electron  beam  we  use  Q  -  2/ n  .  Equation  10  is  written 
intentionally  as  an  integral  over  incident  energy  T  rather  than  local  degraded 
energy  E  in  order  to  facilitate  later  numerical  evaluation  of  Jeo(T). 

The  bremsstrahlung  intensity  viewed  from  the  top  of  the  atmosphere  is  the 
productivity  function  given  in  Eq .  (10)  multiplied  by  an  exponential  attenua¬ 
tion  factor  and  integrated  along  the  appropriate  propagation  path  from  the 
maximum  source  depth  zc  to  the  top  of  the  atmosphere  (see  also  Luhmann,  1977): 

z  seen 

fv(X,n)  ■  /  c  dz  sec  n  P„(2,X)  exp  (~u  z  sec  n)  (11) 

a  o  a 

o 

where  u  *  u  (X)  is  the  mass  absorption  coefficient  which  for  1  keV  <  X  <  30 
a  a 

keV  can  be  approximated 

u  (X)  -  4.28  x  10^  X  cm^/g  (12) 

a 

for  a  nitrogen/oxygen  atmosphere  [  f rom  the  data  of  Hubbell  (1971)]. 

It  is  convenient,  for  our  purposes,  to  evaluate  equations  (10)  and  (11) 
for  discrete  energy  levels  X  *  ;  i  *  1 ,  m.  Also,  we  choose  to  evaluate 
f  (T)  at  these  same  energy  levels  T  »  1-1,  m.  It  is  necessary  at  this 

point  to  specify  the  relationship  between  energy  levels  and  differential 
energy  Intervals  AK^  in  order  to  evaluate  equation  (10)  numerically.  Here  we 
apply  a  geometric  progression  of  energy  levels  which  is  characteristic  of 


current  instrumentation.  Namely, 


where  A  is  typically  between  0.15  and  0.20.  In  doing  this,  tabular  values  of 


4>  ,  computed  for  a  range  of  zenith  angles,  are  specific  to  the  energy 

i  »J 

increments  K,  but  are  otherwise  universal.  Examples  of  values  computed 

J  i  »J 

for  n  *  0°  and  n  *  80°  are  shown  in  Tables  2  and  3  for  values  of  iC j  listed  in 
Table  1.  These  values  are  also  displayed  in  graphic  form  in  Figure  3.  The 
plot  depicts  the  contribution  of  unit  incident  flux  of  electrons  at  energy  T 
to  observable  X-ray  flux  at  energy  K.  Note,  for  example,  the  characteristic 
turnover  of  the  production  function  at  low  X-ray  energies  for  high  incident 
electron  energy.  This  is  due  to  re-absorption  of  low  energy  X-rays  which  are 
emitted  from  deep  within  the  atmosphere  by  deeply  penetrating  high  energy 
electrons.  It  is  useful  to  note  that  a  reasonable  approximation  to  the  pro¬ 
duction  function,  neglecting  soft  X-ray  reabsorption,  is 


i » j 


4.86  x  10 


-6 


(T, 


Kj_)  T. 


k: 

j 


1.5 


(14; 


Of  course,  this  approximate  form  applies  equally  well  to  energy  levels  other 
than  those  specified  in  equation  (13).  A  discussion  of  aspect-angle  effects 
such  as  limb-brightening  or  darkening  appears  in  an  Appendix. 


Table  3.  X-ray  Productivity  for  80°  Zenith  Angle 


<t>  (K,T)/(T- 


Ill .  Maximum  Entropy  Deconvolution 

The  bremsscrahlung  production  function  computed  in  the  previous  section 

(or  any  other  appropriate  kernel  .)  can  be  used  to  compute  the  observable 

i  ij 

X-ray  flux  at  satellite  altitudes  caused  by  incident  electron  flux  with  spec¬ 
trum  fe(T).  In  this  section  a  deconvolution  scheme  is  presented  which  enables 
the  computation  of  the  incident  electron  spectrum  based  on  discrete  observa¬ 
tions  of  the  bremsscrahlung  flux.  The  technique  presented  here  resembles 
similar  maximum-entropy  schemes  which  have  been  developed  for  image  restora¬ 
tion  (see,  for  example,  Gull  and  Daniell,  1978;  Wilczek  and  Drapatz,  1985). 
The  similarity  in  these  data  analysis  techniques  is  to  be  expected,  since 
entropy  maximization  generally  provides  the  smoothest  (i.e.,  most  cautious) 
estimate  of  a  distribution  among  a  set  of  distributions  which  are  allowed 
under  constraints  imposed  by  observations.  Below  we  describe  the  constraint 
equations  appropriate  to  Che  bremsscrahlung  problem  and  develop  a  maximum 
entropy  deconvolution  scheme  for  the  incident  electron  flux  spectrum. 

As  in  Section  II,  we  assume  that  discrete  observations  of  the  bremsstrah- 
lung  X-ray  flux  fx(Xk)  are  available,  and  chat  these  measurements  are  related 
to  the  incident  electron  flux  through  the  following  convolution  function: 


xk 


n 

I 

1-1 


'i.k 


ei 


+  e, 


(15) 


Here  the  e^  terras  represent  errors,  which  we  will  assume  to  arise  from  sta¬ 
tistical  counting  deviations.  This  assumption  requires  an  estimate  of  the 

standard  deviations  of  e^,  which  we  call  o^.  Gull  and  Daniell  (1978)  suggest 

2  2 
the  application  of  a  ^  constraint  to  the  errors  e^  when  their  variances  ok 

are  known.  That  is 
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The  constraint  iT  can  be  chosen  to  reflect  a  confidence  level  p  by  setting 

■y  2 

A  to  the  value  of  the  inverse  x  distribution  function  with  m  degrees  of 
freedom. 

Finally,  the  entropy  of  the  incident  electron  distribution  is  defined  by 
the  Shannon  form. 
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We  wish  to  find  the  distribution  fgi  which  maximizes  equation  (17)  under  the 
constraints  of  equations  (,5)  and  (16).  The  Lagrangian  of  this  system  has  the 
form 
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where  f eg  represents  an  "average"  flux  which  can  be  estimated  from  the  inte¬ 
grated  X-ray  flux,  then  adjusted  in  later  iterations  using  a  continuation 
technique.  Differentiating  the  Lagrangian  with  respect  to  fgi  and  ek,  and 
setting  the  partial  derivatives  equal  to  zero,  we  get  the  following  set  ot  n  + 
rn  equations 
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Solving  for  £ 
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and  substituting  fgi  and  e^  into  equations  (15)  and  (16)  we  get  the  set  of  m  + 
1  nonlinear  equations  for  the  Lagrangian  coefficients  0^ 
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The  system  of  equations  (22)  and  (23)  can  be  solved  iteratively  using  a  non¬ 
linear  Newton's  method.  The  Jacobian  of  the  system  is  given  by  = 

3F,  /  38. ,  where 
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The  iteration  equation  is  then 


7  c#1*  -  6<0) )  -  t  -  0  (25) 

which  represents  a  linear  system  in  8^ 1 ^ .  Following  Wilczek  and  Drapatz 
(1985),  the  starting  points  are  calculated  from  the  initialization  equations: 
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In  test  cases  with  16  discrete  measurement  channels  and  f  .  /  o,  ■  Id  witn  a 

XK.  k. 

confidence  level  of  0.9,  this  scheme  has  been  found  to  converge  in  approxi¬ 
mately  ten  iterations.  Some  examples  are  discussed  in  the  following  sec¬ 
tion.  It  should  be  noted  that  the  solution  set  ot  equation  (25)  allows  an 
easy  "plug-in"  of  different  forms  for  the  bremss t rah  1 ung  production  function 
Thus,  the  deconvolution  scheme  can  be  modified  easily  if  new  or  better 


forms  for  the  production  function  become  available. 


In  this  section,  the  deconvolution  scheme  described  eariier  is  applied  to 
bremsstrahlung  X-ray  spectra  which  are  characteristic  of  the  earth's  aurora. 
The  first  examples  are  "ideal"  in  chat  they  depict  the  ability  of  the  deconvo¬ 
lution  technique  to  reproduce  the  shape  of  Maxwellian  and  Gaussian  electron 
spectra.  These  examples  do  not  include  instrumental  effects,  although  random 
noise  is  applied  to  the  "Ideal"  signal  in  order  to  demonstrate  the  effects  of 
noise  on  the  numerical  technique.  These  first  examples  are  meant  to  give  the 
reader  a  "feel"  for  Che  actual  performance  characteristics  of  the  nunerical 
scheme.  The  second  set  of  examples  simulates  Che  response  of  an  actual  in¬ 
strument  (The  Aerospace  scanning  X-ray  spectrometer  SSB/A  on  DMSP-F6)  to 
ideal  spectra.  The  "noise"  added  in  these  cases  is  characteristic  of  counting 
statistics  variations,  and  consequently  is  energy  dependent.  The  final 
examples  compare  spectral  parameters  inferred  from  actual  bremsstrahlung 
measurements  in  the  aurora  with  parameters  observed  directly  by  a  particle 
detector.  These  final  examples  show  the  performance  of  the  deconvolution 
scheme  in  the  "real  world." 

Figures  4a-c  show  examples  of  the  ability  of  the  numerical  scheme  de¬ 
scribed  in  Section  III  to  reproduce  Ideal  spectral  shapes.  Three  examples  are 
shown:  a  l  keV  Maxwellian,  a  10  keV  Maxwellian,  and  a  5  keV  Maxwellian  with 
an  added  2o  10  keV  Gaussian  peak.  In  these  examples,  the  observable  X-ray 
flux  from  each  input  spectrum  is  calculated  using  equation  (2)  and  the 
bremsstrahlung  production  function  derived  in  Section  II.  A  simple  geometric 
progression  of  16  energy  channels  between  1  and  30  keV  is  assumed.  Random 
noise  is  added  to  this  Ideal  signal  at  the  3,  10  and  30%  level.  The  noisy  X- 
ray  mectrum  Is  then  deconvoluted  by  the  method  described  in  Section  III.  A 


continuous  inferred  electron  spectrum  is  computed  for  energies  between  i  and 
LOO  keV.  The  results  shown  in  Figures  4a-c  are  for  ten  iterations  -  no  sub¬ 
stantial  improvement  was  found  at  higher  iterations. 

It  is  obvious  at  a  glance  that  the  deconvolution  works  particularly  well 
for  the  steep  l  keV  spectrum,  even  at  high  noise  level.  This  characteristic 
is  easy  to  understand  if  one  thinks  of  the  input  spectrum  as  a  simple  signal 
distribution  -  the  1  keV  spectrum  represents  a  sharper  signal  than  the 
smoother,  harder  spectra.  A  second  obvious  characteristic  of  the  deconvolu¬ 
tion  is  its  difficulty  in  reproducing  the  low-energy  portion  of  very  hard 
(e.g.,  10  keV)  spectra  under  even  modest  noise  conditions.  This  problem  is 
due  primarily  to  the  shape  of  the  bremsstrahlung  production  function  itself. 
It  is  a  characteristic  of  the  bremsstrahlung  process  that  energetic  electrons 
produce  much  higher  fluxes  of  low  energy  X-rays  than  do  low  energy  elec¬ 
trons.  The  low  energy  X-ray  signal  is  partially  masked  by  X-rays  resulting 
from  the  high  energy  portion  of  the  incident  electron  spectrum,  and  informa¬ 
tion  is  lost  due  to  decreased  signai-to-noise  ratio  at  low  energies.  This 
effect  is  also  apparent  in  Figure  4c,  where  the  deconvolution  technique  has 
difficulty  reproducing  the  electron  spectral  shape  below  the  Gaussian  peak 
energy.  All  in  all,  the  deconvolution  scheme  produces  good  results  for  noise 
levels  below  about  10X. 

Figures  5a-c  show  the  results  obtained  for  the  same  input  spectra  as  in 
Figures  4a-c,  but  for  a  response  which  simulates  that  of  an  actual  propor¬ 
tional  counter.  (The  Aerospace  Corporation  scanning  X-ray  spectrometer  car- 
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ries  a  proportional  counter  with  geometric  factor  of  .03  cm  ster  and  energy 
channels  as  listed  in  Table  1.  Measurements  of  auroral  X-ray  spectra  are 
obtained  once  per  second  at  an  altitude  of  830  km.)  Naturally,  for  these 
spectra  the  statistical  "noise"  is  energy-channel  dependent,  with  larger 
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relative  errors  occurring  in  the  higher  energy  channels.  However,  because  of 
the  skewed  form  of  the  bremsstrahlung  production  function  these  errors  mani¬ 
fest  themselves  primarily  in  the  low-energy  portion  of  the  inferred  electron 
spectrum.  Again,  soft  spectra  are  better  represented  than  hard  spectra  or 
spectra  with  peaks,  although  reasonable  results  were  obtained  for  all  of  the 
examples  shown  here.  Of  course,  better  results  would  be  obtained  with  instru¬ 
ments  having  larger  geometric  factors,  with  which  statistical  errors  could  be 
decreased. 

Figures  6a,  b  and  7a,  b  show  results  obtained  using  actual  auroral  X-ray 
data  acquired  by  the  Aerospace  scanning  X-ray  spectrometer.  Fortunately  the 
DMSP-F6  satellite,  launched  in  December  1982  into  a  polar,  circular  orbit  at 
840  km  altitude,  carried  both  the  Aerospace  scanning  auroral  X-ray  spectrome¬ 
ter  and  a  precipitating  auroral  electron  detector.  This  instrumentation 
allows  a  good  test  of  the  computational  procedure  described  in  the  previous 
section.  While  the  Aerospace  X-ray  spectrometer  is  a  raster  scanning  instru¬ 
ment,  the  data  used  in  this  discussion  is  taken  from  periods  when  the  instru¬ 
ment  was  in  a  static  nadir-pointing  mode.  This  permits  a  more  detailed  com¬ 
parison  with  the  precipitating  electron  data.  The  electron  data  is  from  the 
J-package  instrument  provided  by  the  Air  Force  Geophysics  Laboratory,  which 
provides  a  precipitating  electron  spectrum  once  per  second  over  the  energy 
range  30  eV  <  T  <  30  keV.  Only  the  data  from  T  >_  2  keV  will  be  used  here. 

In  order  to  compare  the  observed  and  inferred  electron  spectra  we  have 
chosen  to  deal  with  the  fundamental  spectral  parameters;  integral  energy  flux 
(ETOT),  integral  number  flux  (JTOT)  and  characteristic  energy  (EO),  defined 
be  Low : 

30  keV 

ETOT  ■  a  /  dT  J  (T)  T  erg/cm  sec 
2  keV  0 
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Figure  6a.  Plots  of  the  Precipitating  Electron  Energy  Flux  and  Characteristic 
Electron  Energy.  The  solid  line  reveals  spectral  parameters 
computed  from  direct  observations  of  the  precipitating  electrons 
that  cross  indicate  quantities  inferred  from  hreras s t rahl ung  X-ru- 
observations. 
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Figure  bb.  Plots  of  the  Precipitating  Electron  Energy  Flux  for  an  Auroral  Event 
on  29  December  1982.  The  solid  line  reveals  spectral  parameters 
computed  from  direct  observations  of  the  precipitating  electrons  that 
cross  indicate  quantities  inferred  from  bremsstrahlung  X-rav 
observations. 
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Figure  7a.  S3me  as  Figure  5a  but  for  a  Quiet  Aurora  Late  on  29  December  19S 
Note  that  this  example  includes  both  dawn  and  dusk  aurora. 
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Figure  7b.  Sane  as  Figure  5b  but  for  a  Quiet  Aurora  Late  on  29  December  1982. 
Nate  that  this  example  includes  both  dawn  and  dusk  aurora. 
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where  a  •  1.602  x  10~^.  These  spectral  parameters  are  also  important  in  tnat 
they  are  used  often  as  the  basis  of  other  empirically  derived  relations,  such 
as  ionospheric  conductance. 

Two  samples  of  data  were  chosen  for  study.  The  first,  acquired  on  29 
December  1982  between  47320-47460  seconds  UT  in  the  dusk  sector,  is  character¬ 
istic  of  an  auroral  surge.  The  second,  acquired  on  29  December  1982  between 
80200  -  81500  seconds  UT  is  more  characteristic  of  a  quiet  aurora,  and  in¬ 
cludes  auroral  data  from  both  the  dusk  and  dawn  sectors.  The  first  event 
occurred  in  darkness,  while  the  second  was  in  full  sunlight. 

Figures  6a,  b  show  the  observed  and  inferred  electron  spectral  parameters 
from  the  active  auroral  surge  event.  The  satellite  proceeded  from  high  toward 
low  latitude  during  this  event.  The  DMSP  optical  imagery  revealed  that  the 
satellite  passed  over  a  bright,  narrow  arc  at  high  latitudes  before  encounter¬ 
ing  a  broad  core  of  a  westward-travelling  surge  between  47400-4745U  UT.  The 
precipitating  electron  data  (solid  line)  indicates  an  energy  deposition  rate 
of  about  20  erg/cm^  sec  in  the  poleward  arc  and  10-15  erg/cra^  sec  within  the 
core,  with  a  low  level  of  precipitation  separating  the  two.  Observed  charac¬ 
teristic  energies  were  -  15  keV  in  the  arc,  5-6  keV  in  the  gap  and  near  ID  keV 
in  the  core.  The  quantities  inferred  from  the  X-ray  observations  are  plotted 
as  crosses,  with  the  horizontal  bar  indicating  the  effective  spatial  resolu¬ 
tion  of  the  detector  and  the  vertical  bar  depicting  uncertainties  attributable 
to  statistical  counting  errors  in  the  raw  X-ray  spectral  measurement.  The 
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inferred  quantities  are,  by  and  Large,  quite  representative  of  the  real 


world.  The  ability  of  the  X-ray  technique  to  reproduce  the  characteristic 
incident  electron  energy  is  particularly  notable.  The  spatial  distribution 
and  magnitude  of  energy  flux  is  also  well  represented  except  for  a  slight 
underestimate  of  the  energy  flux  within  the  poleward  arc.  This  underestimate 
is  probably  due  to  incomplete  beam  filling  of  the  detector,  pointing  out  the 
future  need  for  well  collimated  instrumentation.  Note  that  in  Figure  6  the 
electron  and  X-ray  data  have  been  offset  in  time  by  15  seconds  in  order  to 
account  for  magnetic  field  line  tilt  between  the  satellite  and  the  atmosphere. 

The  second  example,  shown  in  Figures  7a,  b  depicts  much  quieter  auroral 
conditions,  with  the  energy  deposition  barely  exceeding  1  erg/cm ^  sec.  Note 
that  the  energy  flux  scale  in  Figure  7a  is  one  order  of  magnitude  lower  than 
that  of  Figure  6a.  The  dusk  sector  is  comprised  of  a  number  of  detached  arcs 
with  energy  flux  less  than  1  erg/cm^  sec  and  characteristic  energies  of  about 
3  keV.  The  dawn  sector  is  more  spatially  continuous  with  energy  deposition  of 

*5 

about  1  erg/cm-  sec  at  a  characteristic  energy  of  about  10  keV.  Again,  the 
spectral  quantities  inferred  from  the  X-ray  data  are  in  very  good  agreement 
with  the  observed  values,  particularly  the  characteristic  energies.  The 
poleward  edge  of  the  dawn  aurora  is  underestimated  in  flux,  probably  due  to 
decreasing  characteristic  energy  in  that  region. 

These  two  examples,  which  represent  quite  different  levels  of  auroral 
activity,  show  that  the  X-ray  technique  provides  a  reasonable  representation 
of  the  Important  incident  electron  spectral  parameters,  namely  integral  energy 
flux,  integral  number  flux  and  characteristic  energy.  These  parameters  are 
interesting  in  themselves  since  the  source  of  plasma  in  different  types  of 
auroral  forms  is  currently  a  topic  of  theoretical  interest.  Imaging  X-ray 
detectors  can  provide  large-scale  two-dimensional  maps  of  electron  spectral 
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V.  Summary 

A  numerical  technique  has  been  developed  which  provides  a  statistically 
optimal  estimate  of  incident  auroral  electron  spectra  based  on  satellite 
observations  of  bremsstrahlung  X-ray  spectra.  The  technique  utilizes  a  maxi¬ 
mum-entropy  estimate  for  the  inferred  electron  spectrum  based  on  discrete 
observations  of  the  bremsstrahlung  spectrum.  An  integral  convolution  of  the 
electron  distribution  by  a  bremsstrahlung  production  function  is  assumed,  and 
measurement  errors  are  incorporated  in  the  model.  The  measurement  errors  are 

constrained  in  the  numerical  technique  by  estimates  of  the  variance  due  to 

2 

counting  statistics,  and  a  specified  confidence  level  for  the  x  function.  An 
approximate  form  of  the  bremsstrahlung  production  function  is  computed  in  the 
test;  however,  the  nunerical  deconvolution  scheme  does  not  rely  on  the  partic¬ 
ular  form  derived  in  this  report  .  The  deconvolution  scheme  therefore  can  be 
used  as  a  testbed  for  other  more  accurate  representations  of  bremsstrahlung 
production. 

A  number  of  examples  are  shown  which  characterize  the  performance  of  the  ■ 
numerical  model  both  for  ideal  cases  as  well  as  for  actual  data.  Comparisons 
between  inferred  and  observed  electron  spectra  show  that  while  current  instru¬ 
mentation  suffers  from  relatively  poor  signal-to-noise  ratio,  general  charac¬ 
teristics  of  auroral  electron  precipitation  can  be  measured  remotely  by  satel¬ 
lite-borne  X-ray  spectrometers. 
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Appendix 


The  expression  relating  observable  bremsstrahiung  flux  to  volume  emissiv- 
ity,  shown  below,  indicates  two  angular  effects  of  importance: 


«X,  n)  -  / 


c  sec  n 


dZ  sec  n  P  (x,Z)  exp(-u  z  sec  n) 
3  a 


The  first  effect  is  purely  geometrical,  arising  from  the  f ield-of-view  fore¬ 
shortening  at  increasing  zenith  or  elevation  angle.  This  effect,  along  with 
the  concomitant  integration  over  a  longer  emission  column  driven  by  the  sec  n 
terra,  leads  to  limb  brightening.  The  second  effect  is  increased  attenuation 
due  to  propagation  path  lengthening,  appearing  as  the  exponential  term.  This 

obviously  leads  to  limb  darkening.  Since  the  absorption  coefficient  u  is  a 

a 

strong  function  of  X-ray  energy  one  might  expect  different  trade-offs  between 
limb  brightening  and  limb  darkening  at  different  elevation  angles  and  ener¬ 
gies.  This  is  indeed  the  case. 

Figure  8  shows  a  representative  profile  of  atmospheric  depch  in  altitude 
(Berger  and  Seltzer,  1972)  which  we  can  apply  in  demonstrating  the  effects  of 
limb  brightening.  Figures  9a-d  show  the  ratio  of  the  bremsstrahiung  intensity 
at  the  top  of  the  atmosphere  to  that  at  the  source  altitude  (a  parameter)  as  a 
function  of  zenith  or  elevation  angle  for  several  energies  covering  the 
auroral  range.  For  X-rays  at  20  keV  (Figure  9a),  limb  brightening  is  apparent 
for  all  source  altitudes  from  90-200  km.  The  brightening  is  about  a  factor  of 
five  at  80®  elevation  angle.  Substantial  limb  brightening  is  also  apparent  at 
X  *  6  keV,  although  less  brightening  occurs  for  X-rays  from  the  lowest  source 
altitude.  Lower  energy  X-rays,  at  2  keV  and  1  keV,  show  significant  limb 
darkening  effects  for  low  altitude  sources,  while  maintaining  limb  brightening 
for  higher  source  altitudes.  The  net  effect  of  the  limb  is  an  integral  over 
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Figure  9a  and  b.  Plots  of  the  Ratio  of  Bremsst rahl ung  Intensity  at  the 

Top  of  the  Atmosphere  To  That  at  the  Source  Altitude  as 
a  Function  of  Zenith  Angle  for  (a)  20  keV  and  (b)  6  keV 
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Figure  9c  and  d. 


Plots  of  the  Ritio  of  Bremsstrahlung  Intensity  at  the 
Top  of  the  Atmosphere  To  That  at  the  Source  Altitude  as 
a  Function  of  Zenith  Angle  for  (c)  2  keV  and  (d)  1  keV 


all  source  altitudes,  which  typically  results  in  limb  brightening  at  all 


energies  since  low  energy  X-rays  are  generally  emitted  at  higher  altitudes, 


Note  that  the  effects  of  limb  brightening  or  darkening  do  not  become  signifi¬ 


cant  until  elevation  angles  exceed  about  70®. 


LABORATORY  OPERATIONS 
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The  Aerospace  Corporation  functions  as  an  'architect -engineer "  for 
national  security  projects,  specializing  in  advanced  nllltary  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff’s  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  therraomechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exciraer  laser  development  including  chemical  kinetics, 
spectroscopy ,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications ; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radlometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth’s  atmosphere,  ionosphere  and  magnetospher** ; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 
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